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Abstract 

The notion of brightness is efficiently conveyed in geometric optics as density of rays in phase 
space. Wigner has introduced his famous distribution in quantum mechanics as a quasi-probabihty 
density of a quantum system in phase space. Naturally, the same formalism can be used to 
represent light including all the wave phenomena. It provides a natural framework for radiation 
propagation and optics matching by transferring the familiar 'baggage' of accelerator physics 
function, emittance, phase space transforms, etc.) to synchrotron radiation. This paper details 
many of the properties of the Wigner distribution and provides examples of how its use enables 
physically insightful description of partially coherent synchrotron radiation in phase space. 

PACS numbers: 41.60.Ap 
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I. INTRODUCTION 



The concept of phase space plays an important role in accelerator physics. Useful tools 
such as Twiss parameters, emittance, phase space propagation have been in long use in 
the accelerator community. The extension of classical phase space concept to synchrotron 
radiation is straightforward for geometric optics, applicable for incoherent radiation. The 
Wigner distribution, or Wigner distribution function (WDF), was recognized to be a general 
framework to represent quantum [Ij and therefore wave phenomena in phase space 0[3|. The 
approach allows light characterization of arbitrary degree of coherence [4J and polarization 
[5] in phase space, though its application by the accelerator community has so far been 
mostly limited to simplest cases of Gaussian or Gauss-Schell beams [6^ 7J. This provides 
a set of useful analytical expressions for quick estimates of performance of modern x-ray 
sources with improved coherence properties. In particular, the concepts of the diffraction 
limit and brightness have been extended to cover partially coherent radiation cases following 
the notion of Gaussian distributions in the phase space for both the synchrotron radiation 
of the undulator central cone and density of electrons. More detailed approach to coherent 
or partially coherent sources inevitably calls on physically rigorous wave description of the 
radiation, either through using cross-spectral density [8j or the Wigner distribution. In 
particular, since neither undulator radiation, nor electron distribution in the phase need to 
be Gaussian, the general framework becomes essential to be able to describe the performance 
of x-ray sources with improved coherence. 

The Wigner distribution function provides a natural and elegant description of the prop- 
erties of light, and can serve as a useful tool in accelerator and x-ray beamline design 
including electron to x-ray beam matching, light propagation, fully accounting for arbitrary 
polarization and coherence properties of radiation. The intuitive picture provided by the 
WDF, being the phase space density of light or generalized brightness, is particularly ap- 
pealing to the accelerator community trained to view many aspects of the beam dynamics 
in phase space. Not only can the WDF be readily computed from the first principles, the 
first measurement of x-ray Wigner distribution has been reported in the literature [9j . The 
knowledge of the Wigner distribution represents the entirety of what can be known about the 
radiation and its importance will only increase with advent of more coherent x-ray sources. 

The purpose of this paper is to review many of the useful properties of the Wigner 
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distribution and demonstrate that the WDF can be used with physical insight to describe 
partiaUy coherent synchrotron radiation. In what foUows, the Wigner distribution properties 
are first reviewed in Section [ll] using the language of quantum mechanics. Various examples 
illustrate the physical meaning of the WDF for both pure and mixed quantum states. The 



case of the synchrotron radiation as discussed in Section |TIT] is then viewed as a natural 
extension of the quantum mechanical treatment. Coherence and dispersion properties of light 
as conveniently conveyed by the WDF are emphasized. A special attention is given to light 
polarization, being an important characteristic of synchrotron radiation. Practical matters 
of computing the WDF are covered in Section |IV| which outlines the general procedure 
for obtaining the Wigner distribution first for a single electron, and then extending the 
result to include electron bunches of Energy Recovery Linac as an example. Since neither 
synchrotron radiation nor electron beam in this case have Gaussian phase space density, 
some consideration is given to generalizing the concepts of emittance and brightness to 
describe non-Gaussian distributions. 



II. WIGNER DISTRIBUTION IN QUANTUM MECHANICS 

The Wigner distribution, initially introduced to account for quantum phenomena in sta- 
tistical mechanics [Ij, provides a convenient description of a quantum mechanical system 
in phase space. The Wigner distribution itself does not possess any new information not 
already contained in quantum state itself, which is fully described (together with its com- 
plete time evolution through Hamiltonian %) either by a pure state ij) or more generally for 
a mixed state by its density matrix p = with state weights "^Pj = 1. The 

utility of the Wigner distribution is in convenient and visual representation of the quantum 
system (and by extension wave optics phenomena) in terms of quasi-probability of having 
both phase space quantities (e.g. x^p). Such characterization, being very familiar to ac- 
celerator physicists, is a natural framework of description for a unified phenomena of both 
classical and wave nature reusing many of the concepts from the accelerator field (emittance, 
/3- function, phase space propagation, brightness, etc.). Quasi-probability refers to the fact 
that while the Wigner distribution is normalized to 1 and is used to compute averages of var- 
ious quanitites as expected for a probability density function, the function can take on local 
negative values. This deviation from non-negativity is essential for general quantum or wave 
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phenomena where position and momentum operators do not commute and the uncertainty 
principle must hold contrary to the classical description. Nevertheless, this non-positivity 
does not preclude measurement of Wigner distribution using tomography techniques P . 

The properties of the Wigner distribution have been studied extensively in the context of 
quantum mechanics [1TI[I2], wave optics [HEJIIS] and signal processing p!4l - [T6] . To provide 
a suitable context, the properties of the WDF are reviewed in this section. For simplicity, 
we limit our consideration here to a ID scalar wavefunction, Extension to higher 



dimensions and polarization as required for synchrotron radiation is detailed in Section IV 



A. Pure quantum state 



First, we consider a pure quantum state Particularly insightful definition of the Wigner 
distribution can be given in Dirac notation: 

W{x,p) = J {iP\x + ^){x+ ^\p){p\x - f )(x - f It/;) dx\ (1) 

(The integration here and elsewhere in this paper is taken over the entire range — oc to +oc 
unless stated otherwise.) The integrand is the quantum equivalent of a classical phase-space 
trajectory as seen by reading Dirac brackets from right to left: (1) the probability amplitude 
for a particle in state i/j to have a position {x — (2) the amplitude for a particle with 
position [x — ^) to have momentum (3) the amplitude for a particle with momentum p 
to have position (x + ^)] and finally (4) the amplitude for a particle with position (x + y) 
to (still) be in the state ij). The integration over the entire space x' therefore creates a 
superposition of all possible quantum trajectories of state ip^ which interfere constructively 
and destructively, providing a quasi-probability distribution in phase space [17]. Using a 
well known identity (with h = 27Th the Planck constant) 

^ ipx 

we rewrite 

{x + ^\p){p\x-^) = le"^. 
Eq.[T]then assumes its most frequently quoted form 

W{x,p) = ^J^*{x+^)ij{x-f)e^dx', (2) 
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where + = + y) and ip{x — = (x — I^IV^). 

In the same spirit, Eq.[T]can be rewritten in terms of integration over the entire momentum 
space, 

Wix,p) = J {^\p+P-){p+P-\x){x\p-i^){p-i^\ij)dp\ (3) 

leading to an equivalent definition of the Wigner distribution function now in terms of 
momentum representation of the state ^{p) = {p\ip)' 

W{x,p) = ^ J ^^{p+^-)^{p-^-)e-'^dp'. (4) 

The momentum and position representations, il){x) and if'(p), are related via the Fourier 
transform 

^ 7h "^ipy^' dp. 

A summary of the main properties of the Wigner distribution function is given below. 
Properties that are revisited later for a more general case of a mixed state are denoted by 
an asterisk (*). 

Property 1 (Realness) 

W{x,p) G M. (6) 
This property follows from ^{x^p) = W{x^p). 

Property 2 (Normalization and Marginals*) The WDF is normalized to 1 with its 
projections (or marginals) corresponding to nonnegative probability densities in either po- 
sition or momentum 

r r 

W{x^p) dxdp = 1, 



j Wix,p)dp=\ijix)\^ (7) 
W{x,p) dx = \^{p)\'^. 

The proof is by substitution of the Wigner definition into Eqs. [7]and then using the identify 

J e'fda = hS{b). (8) 
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Property 3 (Boundness) 



m.,p)\ < 2 = ^. (9) 



This property can be proven using the Cauchy-Schwarz inequahty on the definition of the 
Wigner function. 

It is iUustrative to consider when the WDF assumes ±| extrema. An arbitrary wave- 
function %l){x) can be written in terms of even ip(.{—x) = '0e(^) and odd ipo{—x) = — '0o(^) 
parts: = ipd^) + '0o(^)- Then, the WDF at the origin becomes 



2 f 

= — / iIj^{x)iIj{—x) dx. 



This can be written in terms of the wavefunction's even and odd parts: 

W{0,0) = (|^e(a;)r - m^)?) dx. (10) 



As can be seen from Eq. 10 and the wavefunction normahzation, W{0^ 0) = ^ if ip{x) is even 



and W{0^ 0) = — | if is odd and vice versa 

Property 4 (Expectation values) The expectation value of an operator A can be found 
from its phase-space representation function A{x^p) according to 

(i) = J J A{x,p)W{x,p)dxdp, (11) 

where W{x^p) acts as a phase-space probabihty density The function A{x^p) and operator 
A satisfy the foUowing relationships [llj 

A{x,p) = J {x-^\A\x + ^)e'^dx', (12) 
= A{^-^,p)e''^dp. (13) 



The pair of Eqs. 12 and 13 is referred to as the Wigner- Weyl transformation [T3 



Refer to [TT] for proof. A practical significance of this property is that any linear com- 
bination of functions of only position operators or only momentum operators correspond 
to the classical phase-space representation given by Eq. [12] where one replaces x ^ x and 
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p ^ p. In particular, the n-th moments of the distribution are readily obtained using 

(.") = /."»'(x.p)<ix*. 

Quantum mechanical correspondence to the classical correlation expectation of position- 
momentum is more involved due to non-commuting nature of the operators. Indeed, the 
operator xp is not Hermitian, i.e. {xpY = p'^x'^ = px xp since [x,j5] = ih 0. As a 
result, the expectation value {xp) is generally complex and {xp) ^ {px). A solution is to 
write + = {xp) = (px), where the symmetric operator is now Hermitian and its 

corresponding phase space function is found from Eq. 12 to be ^{xp + px) xp. Therefore, 

(.p) = (p.) = //.pinx.P)<ix#. 

The above equations allows us to compute the S-matrix of the quantum phase-space 
distribution familiar to accelerator physicists 

|=.T, (14) 

{px) {p^) / \ -ea 67 / \ -o^ 7 



with the usual meaning of emittance e = Vdet S and Twiss parameters satisfying det T = 1. 
The Heisenberg uncertainty principle can then be written as 

> I (15) 

Property 5 (Time evolution) For a time-independent Hamiltonian ?^ = p^ /2m + V{x)^ 

the time evolution for the Wigner distribution W is governed by 

dW _ _p dW 1 
dt m dx ih 

The proof is straightforward using time-dependant Schrodinger equation. Refer to for 
details. In particular, for a linear force F{x) = Fq — kx with a potential energy V{x) = 
Vq — FqX + (Fq, fc, and Vq are arbitrary constants), h drops out from the Eq. 

we recover the classical Liouville's evolution of the phase-space distribution 
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and 



at m ox op 

This property further illustrates the connection to the classical concept of phase space. In 
particular, classical invariants and transformation rules directly carry over to the quantum 
phase space density in case of no or linear forces. 



Property 6 (State cross-correlation*) Cross-correlation of the wavefunction can be re- 
covered from the WDF of a pure state via a Fourier transform 

^(xi)^*(x2)= [w{^,p)e'^'^dp, 

This property is proven by substituting the WDF definition and using the identity ([s]). 
Note the similarity to Weyl's relationship, the Eq. [T3| 



It should be noted that the cross-correlation function of Eq. [18] is just a density matrix 
of a pure state in either position or momentum basis 

(19) 

This connection of the Wigner distribution to the density operator matrix will continue for 
mixed states as discussed later. 



Property 7 (State recovery*) Property [6] allows to recover the wavefunction from the 
WDF modulo a complex constant 



7/;(x)V^*(0) = / W{^,p)endp, 

% (20) 

Property 8 (Integrated product) For two WDFs corresponding to pure states and x 

W{^^){x,p) = J ^*(^ + - i)e~^dx\ 

W^^){x,p) = ^J x*(x + f )x(^ - 1)6"^ rfx, 

the integrated (overlapping) product is related to scalar state product according to 



If 2 1 2 

W^^){x,p)W^^){x,p)dxdp= - j ^\x)x{x)dx = - • (21) 

The proof of this property again involves a substitution of the WDF definition into Eq. [21} 
the use of identity ((sl) and a change of integration variables. 
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Property 9 (Generalized integrated product) To account for Wigner distribution of 
a superposition of quantum states, we introduce a generalized Wigner distribution (now 
generally complex) using 

W(^u^2){x^P) = T I + |)V^2(^ - \)e^dx'. (22) 



Note that W^*^^ ^ ^(V^2,V^i)- The standard Wigner distribution W(^^) = W(^^^^) is obtained 



by setting '0i(x) = '02(^) = '0(x) in Eq. 22 



Then the integrated overlapping product of two generalized WDFs becomes 

^{^i.H){x^pW{xi.X2){^^P) dx dp = ^ (V^i|x2)(Xi|^2) (23) 
Property 10 (Superposition of states) Consider a superposition of states \iJj) = J2^n 

n 

(either finite or infinite sum). States 0^ need not be orthogonal, but all states are assumed 
normalized: = 1 and ((/>n|0n) = 1- The Wigner distribution can then be written as 

= 5Z 5Z <'^rnW^<Pn.<Pm)' (24) 
n m 

Note that the superposition of states generally leads to appearance of cross-terms VI^((/)^,(/)^) in 
the Wigner distribution. Realness of Wi^^) is readily verified by noting that the oflF-diagonal 
terms in Eq. [24]are complex conjugates of each other (Q^*<^mW^((/)n,(/>m))* = ^m^nW(^^^^^^) = 



oCr^^nW?^ . ^ and therefore their sum must be real. 

For example, consider a stationary Hamiltonian H producing a complete orthogonal basis 
\n) with corresponding energy eigenvalues 

n \n) = En \n) . (25) 

Then, the time evolution of an arbitrary state characterized by the initial vector l^/^o) = 
|'0(t = 0)) adopts the familiar form 

I^W) = E«-l^)^"'^' (26) 

n 

where the expansion coefficients are found in terms of projections of the initial state on the 
eigenbasis = (nlipo) and must satisfy normalization requirement J2n I^^P ^ ^• 
The time evolution of the WDF for |'0(t)) is given by 

W^^){x,p]t) = ^ ^ <a^e W^n,m) {x,p). (27) 



Eq. 27 can be also rewritten using (23) in terms of integrated products of the Wigner 
functions 



n m 

Note that the decomposition of the Wigner distribution for a pure state into an othronormal 
basis generaUy requires the presence of non- vanishing interference cross-terms. We shaU 
revisit this subject when considering incoherent addition of states. 

Property 11 (Gaussian state) A positive WDF can only be realized for a wavefunction 
of the form 

iP{x) = e-(«-'+^-+^), ^{a} > 0, (29) 

leading to a joint Gaussian WDF in position and momentum [12]. We also note in passing 
a well-known fact that a Gaussian state yields the phase-space probability density with the 
smallest rms spread (quantum emittance) e = h/2. 

Property 12 (Convolution) The distribution function obtained by convolving two WDFs 
each corresponding to an arbitrary pure state is everywhere positive. 

This property is introduced in [T9] and is mentioned here for completeness. It should be 
noted that the distribution function so obtained is no longer the quasi-probability suitable 
for finding the expectation values of the original state. For further implications of this 
property, including it being a possible mathematical tool of the concept of measurement in 
quantum mechanics, the reader is directed to the discussion in [19]. 

1. Example: wave packet time evolution in ID potential 

Next, we consider several examples that illustrate the concept of the Wigner distribution. 
The first example shows the phase-space motion of an electron in ID potential depicted in 
Fig. [T] The potential consists of a perfectly refiecting barrier on the left and a simple 
harmonic oscillator (SHO) potential on the right. In units of eV the potential is given by 
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W(^^){x',p')W(m,n){x',P') dx'dp' 



i{En-Em)t 



(28) 




0.5 1 1.5 2 2.5 3 3.5 4 

X (nm) 



FIG. 1. Potential for an electron V{x). The initial state Tp{x,t = 0) and the first 40 energy 
eigenstates tpnix) offset by their eigenvalues are shown. 



(the position x is in nm): 



V{x) = < 



oo, a; < 

0, < X < 2 

10{x - 2)2, x>2 



(30) 



The initial wave packet is described by a Gaussian ip{x^t = 0) (x e ^^^^ with xq = 3 
and Gx = 0.3. Fig. [l] shows ilj{x^t = 0) and first 40 energy eigenstates '0n(^)- The intial 



quantum state is then evolved according to Eq. [26| The motion in the phase space of a 
classical particle with an initial position x = 3 and a zero velocity is shown in Fig. [2] along 
with expectation values for position Xav = {x) and velocity Vav = {p)/m for the quantum 
case. Points on the plot correspond to T/16 time steps where T is the natural period of 
motion of the system, equal to about 3.8 fs in this case. The apparent damping in the 
quantum case corresponds to delocalization of the initial Gaussian wave packet. 

Fig. [3] shows the Wigner distribution of the quantum motion along with its classical 
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01234 01234 
x^i (nm) x^^ (nm) 

FIG. 2. Motion of a classical (left) and quantum (right) electron in phase space. In both cases the 
initial position is 3nm with no velocity. Points show time steps from to T with steps of T/IG, 
where T — 3.8 fs is the natural period of the motion. 

counterpart (marked by x) at the same times as depicted on Fig.[2| The Wigner projections 
K^l?/;)!^ and KpI?/^)!^ are also shown. In addition, the rms emittance of the quantum phase 
space distribution is shown. As can be seen from Fig. [3| the emittance of the initial wave 
packet is increasing when the wave packet reaches the hard reflective potential boundary 
(and more generally when discontinuities in the potential are encountered). Here and in all 
subsequent plots of the Wigner distribution we use the same color map: blue and red colors 
correspond to negative and positive values respectively and the white corresponds to zero. 
It should be noted that despite increase in the phase space area, the mode obviously remains 
pure at all times. 



2. Example: eigenstates of a simple harmonic oscillator 

Another example we consider is the WDF of energy eigenstates for a SHO [20], which 
adopts the same math as Hermite-Gaussian modes in optics (at a waist and a single coor- 
dinate). The well-known eigenstates of the Hamiltonian = + ^re of the form 
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-2/h 



+2/h 




FIG. 3. Time evolution of the Wigner distribution. Solid curves show projections of the WDF, 
probability densities in position and momentum (velocity) respectively. The electron's classical 
counterpart is depicted by x . 



_ 1 

ijn{x) = ^^^e-^H^{x), with n = 0, 1, 2, . . . . (31) 

Hn are Hermite polynomials. For this example we set the usual constants h^m^u) ^ 1 to 
minimize clutter in the expressions. The Wigner distribution is then a radial function in 
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FIG. 4. First three states (left) of a simple harmonic oscillator along with the radial Wigner 
distribution (right). 



W = ^ ^e-" Ln{2r'). (32) 



phase space with = + 

Here are Laguerre polynomials. The first three states along with W{r) are shown in 
Fig. [4} We make a couple of observations regarding Eq. [32 : 



The WDF is the maximum possible value at the origin for even n and the minimum 



possible for odd n in accordance with Eq. 10: VI^(O) = (— l)^/7r or in the regular units 



iy(x = 0,j9 = 0) = (33) 

• The uncertainty in position and momentum or emittance for mode n is 2n + 1 times 
1/2 in the natural units or h/2 in regular units. In optics, e/{h/2) quantity is known 
as the M^-parameter, i.e. 

= 2n + 1. (34) 

In natural units used in this example we have 

(x^) = (/) = e=— . (35) 

Thus, the ground Gaussian state has the smallest possible uncertainty of 1/2 (or /^/2), 
whereas each subsequent excitation adds an additional node to the wavefunction and the 
radial Wigner distribution, and increases the emittance by 1 (or h). 

14 



B. Mixed quantum state 



Generally real quantum systems cannot be described as a superposition of pure modes 
(which is itself a pure mode) instead adopting a statistical language to describe an incoherent 
mixture of pure states characterized via the density operator 

P = '^Pj\^j){^j\^ (36) 
j 

with state probabilities pj adding up to one. When only one coefficient pj = 1 for some j is 
present the formalism is reduced to that of a pure state. Recall that an expectation value 
of an operator A is given in terms of a trace 

{A) = j2pjmm) = MpA). (37) 

j 

Other standard properties are 

Tr(/3) = 1, (38) 
Tr(/32) < 1, (39) 

where the equal sign in Eq. [39] is for a pure state case and less than 1 otherwise. 
The definitions Eq. [2] and [4] are now replaced with 



(40) 



In other words, the WDF of a mixed state is a weighted sum of the WDFs corresponding to 
individual pure states of the density matrix 

Wix,p) = J2PjWi^.)i^,p). (41) 

j 

Since, the Wigner distribution is a quadratic (intensity-like) function of the state, addition 
of individual WDFs corresponds to an incoherent addition, which is to be contrasted with 



coherent superposition of Eq. [24} This language of a pure state vs. a mixed state, incoherent 
addition vs. coherent superposition carries over directly to optics and allows characterization 
of partially coherent sources. 
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Next, we present a new property that relates the WDF to modal purity. We then revisit 
some of the properties introduced earlier to extend them for the mixed state case. Most of 
the properties discussed previously, namely Properties [T] through [7] directly carry over to the 
mixed state case after the necessary modifications to Properties [2} [6j and [7| 

Property 13 (Measure of modal purity) The integrated WDF squared is a measure of 
modal purity 

Tr(p2) = hll W\x,p)dxdp < 1. (42) 



The equal sign in Eq. [42] is for a pure state. 

The proof of this property follows directly from Eq. [41] and Property [8| 
Property 14 (Marginals — Property [2] revisited) 

J W{x,p)dp= {x\p\x) = ^pj\^j{x)\^, 

(43) 

W{x,p)dx= {p\p\p) = ^pj\^j{p)\^. 

3 

Again, this property reinforces the notion of simply adding the intensities (here probability 
densities) for mixed states. 

Property 15 (Density matrix — Property [6] revisited) The density matrix is related 
to the Wigner distribution via a Fourier transform 



p(xi,X2) = (xi|p|x2) = y* W[^^^,p) '^''(ip, 

W{x, e ^ dx, 

which are simply the inverse of the Wigner function definitions 

W{x,p) = j p(a; — y, x + y)e^da;', 



(44) 



(45) 



Property 16 (Mode decomposition — Property [T] revisited) Property [T] to invert 
the wavefunction from its WDF is only applicable for a pure state returning a meaningless 
''wavefunction" otherwise. Since the density matrix is a posit ive-semidefinite Hermitian 
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operator with unit trace, it has an orthonormal basis of eigenstates (j)n whose corresponding 
real eigenvalues > and = 1 : 



TV 



p = Y,K\(t)n){K\. A, >0, 5^A, = 1. (46) 



n=l 



Therefore, from Eq. [45| the WDF can also be written as an incoherent sum of orthogonal 
pure modes 



N 



W-(a;,p) = 5]A„W^(^„)(x,p). (47) 



n=l 



The knowledge of the Wigner function or the corresponding density matrix p allows 
finding the modes and their weights via the standard eigenvector and eigenvalue problem. 



i.e. as seen from multiplying both sides of Eq. 46 by \(j)m) and using orthonormality condition 

p|0m) = A^|(/)^). (48) 

Any convenient complete orthonormal basis can be chosen to represent the density matrix 
p and its eigenstates 0^ [22j. 

Some additional comments about this property are in order. The decomposition given 
by Eq. [47|is distinct from Eq. [41] in that the decomposition yields orthogonal states, which 
is generally not the case in how the mixed state has been originally set up. As a result, the 



modes of the decomposition, Eq. 46, may bear little resemblance to the original preparation 



states of the density matrix (i.e. different mixtures may correspond to the same density 



operator). For the case of a mixed state with orthogonal preparation states, Eq. 41, the 
decomposition recovers the modes and their weights exactly. Small negative eigenvalues 
A^ usually indicate an experimental error in arriving at the density matrix (or the Wigner 
distribution) [22] and can serve as a diagnostics and a self-consistency check. Finally, as is 
generally the case for pure states, the modes 0^ need not be simple in the sense that they 
don't necessarily have a small momentum-position uncertainty (i.e. can be ^ 1). 
As an example, consider a relatively complicated mode of Fig. |3] after one or more oscillation 
periods, which, despite having a large dispersion, is still a pure (fully coherent) mode with 
Tr(p^) = 1. The use of Property 16 recovers just that mode, which itself may have a very 



rich spectrum in some other basis. 
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1. Example: superposition of two states 



In this example we demonstrate the difference between a coherent superposition and 
an incoherent mixture of two Gaussian states. A numerical example of a mixed mode 



decomposition further demonstrates the use of Property 16 



First, let us consider another property useful in this example. 

Property 17 (Quadratic phase) The effect of multiplying an arbitrary wavefunction 
iIj{x) by a pure phase factor of the form g^(co+cix+2C2a: ) arbitrary real coefficients Cq, Ci, 
and C2 

= V^(x)e^(^°+^^"+i^^"'\ (49) 
leads to the WDF related to the original W(^^) via the linear momentum transformation 

W^^^{x,p) = W^^){x,p), with p = p- h{ci + C2x). (50) 

In other words, multiplying the wavefunction by a linear phase factor amounts to a shift 
in momentum, whereas the quadratic phase shift adds a linear correlation (chirp) to the 
momentum vs. position. 

Fig. [5] shows a superposition and an incoherent addition of two Gaussian wave packets. 
In this example, the position is in nm and velocity is in nm/fs (an electron is assumed). 
Coherent superposition is for two Gaussian wave packets with equal weights (the Gaussians 
are nearly orthogonal as seen from the fact that their Wigner distributions don't overlap 
and the Property Isl) . 



Consider a superposition of two states oc + \ip2)- Using Property 10, the Wigner 
distribution has 3 terms 

W^^) oc W^^,) + W^^,) + ty„ (51) 

where the interference term Wi = W(^^^^^^) + W^^^ The interference term is responsible 
for oscillations seen in Fig. [5| and it is easy to show that Wi for two orthogonal states carries 
no energy 

J J W,{x,p)dxdp = 0, if (7/;i|V^2) =0. (52) 

On the other hand, a mixed state p oc \ipi){ipi\ + |'02)('02| has only 2 terms from each 
individual states in its Wigner distribution 

W^p)(xW^^,) + W^^,), (53) 
18 




FIG. 5. The Wigner distribution for coherent (left) and incoherent (right) superposition of 
two quantum states with equal weights. Solid curves show projections of the WDF, prob- 
ability densities in position and velocity respectively. The states are of the form t(;(x) oc 
g-Ax /2cr^^i^[vo x^2dx ^ \ whcrc Ax = X — xq. The parameters {x^^ax^v^^dv / dx} are equal 
to {3nm, 0.3 nm, Onm/fs, 4fs~^} and {2nm, 0.5 nm, 3nm/fs, 1.5 fs~^} for the two Gaussian states 
depicted. 



without the interference term. 



Fig. |6] demonstrates orthogonal mode decomposition using Property [16] for a mixed state 
with two Gaussians of equal weights. As seen from the Fig. [6| recovered modes generally 
reflect the shape of the mixed state and can have large dispersion (emittance) if the mixed 
state itself has had a large dispersion. 
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FIG. 6. The Wigner distribution for a mixed state of two Gaussians of equal weight (left column) 
and corresponding orthogonal mode decomposition (middle and right columns). The prepara- 
tion states are of the form t/j{x) oc e~ '^^xQ^2hdx ^ , where Ax = x — xq. The parameters 
{x{)^(Jx^dv / dx} are equal to {3nm, 0.3 nm, 4fs~^} for fixed and {1-2-3 nm, 0.4 nm, lfs~^} for 
displaced Gaussians. The reconstructed WDF is recovered through = ^^^{(t)i) + -^21^(02) 



2. Example: Gauss-Schell model 

Now we consider a quantum mechanical analog of what is known as a Guass-Schell model 
in optics j23]. Using same natural units of SHO with h^m^uj 1, we rewrite the Wigner 
distribution similar to Eq. [32] in a generalized Gaussian form 

where as previously = Setting 1 recovers a pure Gaussian ground state, 

whereas > 1 corresponds to a mixed state. As previously, Eq. 35 applies for our choice 
of units 
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Next, we use Eq. [44] to recover the density matrix 

p(xi,X2) = _^e-i[^'(^^-^^)'+^"'(^^+^^)'l. (55) 
A more common form of presenting the density matrix is as a ScheU-model source 

X2) = ^/I{xl)^/I{x2)^i{xl - X2), (56) 

where the probabihty density 

1 M 

lix) = pix, x) = —= — e 2^^^ , with a^. = (57) 

and the degree of spatial coherence 

//(xi - X2) = e , with = -j===. (58) 

The function /x(Ax) is bound < |//(Ax)| < 1 with 1 or corresponding to a perfect or 
no phase correlation respectively, is known as a coherence length in optics. E.g. 1 
yields ^ oc (a perfect phase correlation or pure state) whereas 00 gives ^ 

(no phase correlation in the state). 



Decomposition eigen problem (48) can be rewritten as 

{x\p\x'){x'\(l)m) dx' = \m , 

p(x, x')(j)m{x') dx' = A^0^(x). 



/ 



This Fredholm integral equation yields the following spectrum of eigenvalues and eigenfunc- 



tions for the density matrix of Eq. 55 [231 El] 



^m{x) = ^ e "2 Hm{x), (59) 

Am = Aog"^, with m = 0, 1, 2, . . . , (60) 

where 

- 1 
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Thus, the Gauss-ScheU model adopts a particularly simple mode decomposition, which are 
the pure states of a simple harmonic oscillator. Additionally, it can be checked that 



oo 



Tr(p) = J]A^ = 1, (63) 

oo ^ 

T^<P') = E^'m = jfr (64) 

m=0 

Eq. 64 can be a source of confusion in that one may be tempted to equate (phase- 
space area, or emittance, or dispersion) directly to spectral purity Tr(p^) for an arbitrary 
mixed state. This temptation should be resisted since quantum-mechanically the two con- 
cepts, the phase-space uncertainty and the mode purity Tr(p^), are distinct as argued 
previously. To the extent that the Gauss-ScheU model is applicable to a quantum or optical 
system, such blurred interpretation of simultaneously being a measure of dispersion and 
coherence may be justified. A notable exception is when 1, which corresponds to 

both perfect coherence and minimum uncertainty of a pure Gaussian state. We shall see 
later, however, that the synchrotron radiation from an undulator source by a single electron 
is far from a Gaussian and, therefore, such dual interpretation of needs to be rejected 
for the diffraction-limited electron beams. Similarly, the use of the Gauss-ScheU model on a 
distinctly non-Gaussian phase space distribution function has little merit. 



III. WIGNER DISTRIBUTION FOR SYNCHROTRON RADIATION 

The connection of the Wigner distribution to describing partially coherent sources is 
usually made through the cross-spectral density function r(ri,r2,CL;) [4j 

r(ri, r2; c^) = {E{t^,uj)E%T2; u)) . (65) 

Here E{ri^uj) is frequency representation of the electric field, which assumed for now to 
be a scalar function (e.g. linearly polarized light) of 2D transverse coordinates r = (x, y) 
(e.g. the detector plane) and (. . .) means ensemble average (e.g. over electron bunches for 



synchrotron radiation). For Eq. 65 to fully describe coherence properties, the source needs 
to be stationary in that all ensemble averages do not vary with respect to time (or at 
least first and second moments are time-independent, which is a requirement for wide- 
sense stationary processes). The synchrotron radiation with its pulsed bunch structure is 
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generally non-stationary. However, as argued in [8], one can use the cross-spectral density 



in the form of Eq. [65] if individual synchrotron radiation pulses last much longer than their 
coherence time (the time scale of short-term field fiuctuations, inversely related to the source 
bandwidth) , or at ^ NuUJq for an undulator source with A^^ undulator periods and resonant 
(radiation) frequency cjq and electron bunches of at duration. This condition is usually well 
satisfied (though an extension of the formalism can be straightforwardly made to describe 
nearly transform-limited sources in time- frequency domains) . One also typically defines the 
spectral degree of coherence 



/x(ri,r2;c^) = (66) 
Vr(ri,ri;cj)^r(r2,r2;cj) 

The modulus of the spectral degree of coherence ranges from to 1 for incoherent to fully 

coherent sources, < < 1. For a fully coherent radiation, = 1 everywhere. This 

quantity is directly related to the fringe visibility in interference experiments. 

A fully equivalent characterization can of course be made in time domain [26]. In what 
follows, we restrict our treatment to frequency domain, being a more natural choice for 
x-rays. Therefore, the frequency dependence for the functions will be understood while the 
symbol itself will usually be omitted from the expressions, e.g. E{r) = E{r]uj). 

The Wigner distribution for optics is then given by ^ 

W{r,e)=(l] /r(r-^,r+^)e^^^'^rfV, (67) 



2 



^ , , r{d-^,e + ^)e-'''-^ d'e', (68) 

A 

where transverse position r and angle Q = {9x^ 9y) form a conjugate pair similar to position- 
momentum in quantum mechanics (small angle approximation is used throughout). Cross- 
spectral density in position and angular representations are defined according to 

r(ri,r2) = (E(ri)E*(r2)), (69) 

r(ei,e2) = (£:(ei)r(e2)), (70) 

where the angular representation £^(9) of radiation (far field) is related to its spatial repre- 
sentation E{y) via the Fourier transform pair 

£{Q) = \ f E{r)e-'^'-^(fY, 

(71) 

^(r) = i y 8{Q)e'^'-^d?Q. 
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The radiation wavenumber k above is given by = 27r/A = cj/c in terms of wavelength A, 
frequency cj, and the speed of hght c. 

The connection to quantum mechanics now becomes obvious. Refer to Table [1} All the 
properties introduced in the previous section have their counterparts in optics. In particular, 
we note that the radiation wavelength A in wave optics plays a role of Planck constant h 
in quantum mechanics. E.g. geometric optics is recovered in the limit A ^ just as the 
classical behavior can be obtained through h ^ 0. 



TABLE I. Correspondence between quantum and optical formalisms. 



Quantum mechanics 


Wave optics 


State 

Planck constant 
Uncertainty principle 
Conjugate pair 
Density matrix 
Classical mechanics 


h 

e > h/47v 
X < > p 


E{r) 
A 

^x,y > A/47r 

r < > Q 

r(ri,r2) 
A^O 


Field 
Wavelength 
Diffraction limit 
Conjugate pair 
Cross-spectral density 
Geometric optics 


Phase space density 
Normalized to 
Measure of state purity 


W{x,p) 
1 

hJJW^dxdp 


w{r,e) 

spectral flux 

.2 JJW^d^rd^e 
^ (JJWd^rd'^e)'^ 


Spectral brightness 
Normalized to 
Measure of coherence 



Also, the overall degree of coherence //^, which is directly equivalent to Tr(p^) of density 
matrix p in quantum mechanics, can be expressed in terms of the Wigner distribution 
function 

{nW{r,d)d^rd^dy' ^ ^ 

where the denominator, the total flux squared, plays a normalization role so that < /x^ < 1. 

A. Polarized light 

Treatment of polarized light [27] is directly analogous to WDF of spin-| quantum particle 
[28] . which require a 2-component spinor to characterize a state. The reason that a spin-1 
particle (photon) can be described by a 2-component spinor (as opposed to 3) is well known 
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in that only ±h spin projections along the direction of propagation (helicity) are realized 
for a massless particle. 

The Wigner distribution now becomes a 2 x 2 matrix W(r, 9) (complex for off-diagonal 
elements) with components defined according to 



Wkiir, 9) = Q) I {Eliv + ^-^)Ei{v - ^)) e^'^'-^Sv', with I 



x,y. 



(73) 



Generalizing the formahsm of polarized light [271129] . W(r, 9) can be represented as a scalar 
function on the Poincare sphere using 



W{r,e,n) = 0-S(r,9), 
where generalized Stokes parameters are found from 

^^(r, e) = Tv[ajW{r, 9)], with j = 0, 1, 2, 3. 



(74) 



(75) 



(76) 



Here aj are 2x2 Pauli matrices with ctq being an identity matrix, and is a vector mapping 
Stokes parameters onto the Poincare sphere with polar x ^^id azimuthal (j) angles 

/ 1 \ 

sin X cos (f) 
sin X sin (j) 
\ V3C08X I 

The generalized Stokes parameters, which now play a role of a 4-component phase space 
distribution, can we written explicitly in terms of the WDF components 

5'o(r,9) = V^;,(r,9) + W,,(r,9) 
5i(r,9) = W,,(r,9) + l^,,(r,9) 
S2{v,Q) = i[W,y{v,Q)-Wy,{v,Q)] 

In what follows we may occasionally refer to the generalized Stokes parameters as simply 
Wigner distribution functions. These four functions completely characterize radiation of 
arbitrary degree of coherence and polarization in phase space. They can be propagated in 
linear optics just like the scalar WDF. The usual Stokes parameters are found from those 
in Eqs. [TTlby integrating away angles 



(77) 



Sj{r,Q)d^Q, with j = 0,1,2,3. 



(78) 
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Thus, the generahzed Stokes parameters have their usual meaning for polarized light (the 
exact ordering of the Stokes components 1 through 3 may differ in literature): ASo(r, 9) 
represents total intensity in phase space, 5'i(r, 9) represents +45°/ — 45° linearly polarized 
light (for +/— respectively), S'2(r, 9) corresponds to right/left-hand circular polarization, 
and 5^3 (r, 9) to x/^-linear polarization. We note that exact signs here apply only to in- 
tensity projections Sj{r) since the WDF (Sq) is allowed to take on negative values while 
its projections (marginals) are guaranteed to be positive. For example, a Gaussian mode 
with x-polarization will have 5'o(r, 9) = AS3(r, 9) > with other Stokes parameters being 
0, or for left-hand circular polarization 5^2 (r, 9) = — S'o(r, 9). Whereas fully polarized light 
satisfies Sq = ^/sl + + si and partial polarization manifests itself as Sq > V^^f^T^f^T^, 
the generalized Stokes parameter Sq can take on local negative values and deviate from these 
expressions. 

As shown in [27] , the overall degree of coherence for vectorial waves can be written as 



27rA 



2 JJd^rd^e J^^d^QW^{r,e,Q) 



where d'^Q = sinx dxdcj). Or equivalently in terms of generalized Stokes parameters 



(79) 



[JJ So{r,e)d^rd''eY'' 
and explicitly in terms of the WDF components as 



(80) 



A 



2 1 1 (Wl + 2W,yWy, + W^^) d'r d'e 

[JJ{W,:. + Wyy)d^Vd^e]^ 



(81) 



B. Wigner distribution projections 



One of practical limitations of the Wigner distribution is that generally one needs to 
employ four-dimensional arrays as a function of light frequency (and possibly time if the 
temporal structure inside an individual synchrotron pulse is important) times 4 for Stokes 
parameters to represent the radiation fully. In addition to large memory requirements, one 
typically prefers to visualize two-dimensional projections rather than the entire phase space, 
much as it is done in accelerator physics for particle tracking. Here, we mention some of 
the properties of such projected WDFs, limiting our discussion to linearly polarized light for 
simplicity. Important 2D projections of the Wigner distribution are intensity /(x, y), the far 
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field (angular) intensity X{6x^ Oy), x — 6^ and y — Oy phase space projections, Bx{x, 6x) and 

i{x, y) = jj ^{^^ y, ^x, Oy) de, dOy, (82) 

i{e,,ey) = jj w{x,y,e,,ey)dxdy, (83) 

13,{x,e,) = jj W{x,y,e,,ey)dydey, (84) 

Byiy,ey) = jj Wix,y,e,,Oy)dxde,. (85) 

Carrying out the integration and using the identity / e''"'^da = 27rS{b)/k, one obtains 

I{x,y) = {E'^{x,y)E{x,y)), (86) 

B.{x, 0.) = jff r,(x - I, X + |)e^^-^^-rfx^ (87) 



1 

A 

x[d,,dy) = ^E\Q,,dy)E[d,,dy)), (89) 

where r^(xi,X2) = / {E{xi,y)E''{x2,y)) dy and Vy{yi,y2) = / {E{x,yi)E%x,y2)) dx. 

If the radiation modes are separable, i.e. can be written in the form E{x^ y) = (/)x{x)(l)y{y) 
(for example Hermite-Gaussian modes), then all the properties discussed in Section [lT| for 
two-dimensional WDF in quantum mechanics apply to the Wigner 2D projections after 
normalization Wx^y = Bx^y/J^^ where = JJ I{x^y) dx dy is the total (spectral) flux. It 
includes the interpretation of the A JJ W^(x, Ox) dxdOx and a similar expression for y-plane 
to be a measure of coherence IJi^x^y (the analog of the Tr(p^) in quantum mechanics). On 
the other hand, for non-separable radiation flelds (e.g. general radially symmetric modes), 
the same interpretation of A JJ W^{x^ 9x) dxdOx = f^^x cannot be made. 

Nevertheless, for simple linear optics without coupling of x, ^-planes (drifts and lenses), 
the WDF projections can be propagated in the same way as the full four-dimensional Wigner 
distribution. We also note that a pure mode with symmetric flelds E{—x^—y) = E{x^y)^ 
which are of a practical importance to synchrotron radiation, the on-axis 2D brightness takes 
on the possible maximum value 

max{i3,,^(0,0)} = ^J*, (90) 
where is the total spectral flux contained in the mode. 
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C. Light propagation 



One of the strong appeals of the Wigner distribution function is in its natural propagation 
for linear optics, which is entirely similar to the classical phase space evlolution. As a result, 
the formalism developed in the accelerator physics for classical phase space distributions 
can be directly carried over to (partially) coherent synchrotron radiation. In analogy to 
Property [5| the local values of WDF stay constant on phase space trajectories subject to 
classical transformation in drifts and lenses along the longitudinal position z 



iy(r(^2),e(^2)) = H/(r(^i),e(^i)), 



(91) 



where 



( x\ 

y 



M(zi ^ Z2) 



z2 



y 



with detM 



(92) 



Similarly, for a decoupled in x, y-plane transport, the 2D projections of the WDF follow the 
classical transformation 



By{y{z2),6y{z2)) = By{y{zi),6y{yi)). 



(93) 
(94) 



Drift or lens transformations lead to a rotated or sheered WDF, and since the projections 
of the WDF are accessible for measurement, this allows a reconstruction of the Wigner distri- 
bution through tomography, similar to the use of tomography in phase space reconstruction 
in accelerators. 

The introduction of spatial filters (e.g. a pinhole or a slit) naturally leads to diflFraction 
phenomena and the Wigner distribution gets altered in a non-trivial way. Whereas, the 
electric field after an aperture with transmission t(r) is simply E{r) £'(r)t(r), the Wigner 
distribution is given by the convolution of the angular variables 9 of the input Wigner 
function with that of the spatial filter [5] 

W{r, 9) ^ / W{r, e)Wt{r, 9 - 9^ d^Q', (95) 



where 
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IV. NUMERICAL EVALUATION 



In this section we discuss practical matters pertaining to computing the Wigner distribu- 
tion for undulator radiation. As we shaU see, as long as the effects of cosh-dependence of the 
undulator fields can be ignored, the synchrotron radiation in phase space can be obtained 
by a convolution of the WDF from a single electron with that of the entire electron beam 
phase space. This is a consequence of the well-known fact that the electrons in a single 
bunch do not interfere with each other unless there is a microbunching structure on the 
wavelength scale. In which case, the computation of radiation fields proceeds diflFerently. In 
what follows, we limit our examples to the "electron only interferes with itself" scenario, as 
applicable for non- free-electron-laser (non-FEL) emission regimes. 



A. Radiation field generation 

Calculation of radiation fields is well established, e.g. see [30l[3l]. The frequency repre- 
sentation of the electric field is given by 



E r;cj = / - 

^ ' ^ 47reoC J R 



e'^'^^^^'^'^T, (96) 



for an observer at r, the position vector from the observer to the electron R = r — Te, with 
re(r) being the electron's trajectory as a function of time r, the velocity |3 = c~^dr^/dT^ 



and the unit vector n = R/i? with R = |R|. The expression 96 is exact and is convenient 
for numerical evaluation in that once the trajectory re(r) is found, the integral evaluation is 
direct. We assume transversality of the field, i.e. E ^ [Ex-, Ey^ 0). An expression with parax- 
ial approximation for the field can be obtained however, it represents little advantage 
over the exact expression for numerical work. A simulation tool has been developed that 
solves for the electron trajectory in arbitrary field configuration and evaluates the radiation 
integral, Eq. [96l 

The undulator magnetic fields are taken of the usual form 

Bx = Box sin(/c^z) cosh(/c^x). 

By = Boy co8{kuz) cosh(/c^2/), (97) 
Bz = Box co8{kuz) 8mh{kux) - B^y sin(/c^z) sinh(/c^?/). 
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Here A^^ = 2Ti/ku is the undulator period, B^^^y are the maximum magnetic fields in both 
planes with B^^ = for a conventionally oriented planar undulator. The total undulator 
length is taken to be = NuX^ and the relation to the undulator K parameter is via the 
usual Kx^y = eSox,2/A^i/27rmeC. To ensure on- axis orbits with no net defiection, the undulator 
fields are 1/4 and 3/4 of their nominal values for the first and second period halves on either 
undulator end. 

For numerical evaluation of the Wigner distribution function, the Fourier transform of 



Eqs. [87| [88]is replaced with its discrete analog. A detector is placed at an arbitrary position 
z downstream of the undulator, and the electric field is evaluated on a transverse grid of 
positions r^^^^ = (xk^yi^z). The phase space distribution is then typically back-propagated 
to the undulator center using the usual transforms. It should be noted that the discrete 
Fourier transform can suffer from aliasing problems and, in order to avoid this problem, the 
maximum angular extent of the radiation must be within n/kAx^y^ where Ax^y is the grid 
size of the radiation field sampling. To avoid very small grid sizes, it is convenient to use 
Property [17] to first remove the quadratic phase present in the radiation pattern. This is 
equivalent to introduction of a perfect thin lens, which is subsequently removed after the 
WDF is evaluated but prior to phase-space propagation to a point of interest. 



B. Electron bunch effect 

The effect of adding radiation from many electrons is equivalent to an earlier considered 
example of superposition from two quantum states. For any two electrons in the bunch, the 
electric field will differ by a phase factor e^^^^ , where tj represents time of the electron inside 
the bunch. It is easy to see that the interference term of Eq. [51] averages out to since it 
contains essentially random phase factors e^^^^^^"^'^) inside the averaging brackets. In other 
words, the uncertain phase relationship between the two electrons on the optical scale leads 
to a density matrix case analogue where the interference term drops out and the WDF is 
simply an incoherent sum over all the electrons. 

Therefore, if the Wigner distribution function of a single electron does not change its 
shape, but simply shifts in position and angle, as one would expect for an undulator in 
which the trajectories remains linear with small position and angle offsets, the overall radi- 
ation pattern is just a convolution (summation) of the electron distribution in phase space 
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with the WDF of a single electron. Additional effects arise for either segmented undulator 
(with focusing between the segments), due to the vertical focusing of a planar (horizon- 
tally deflecting) undulator, or due to the effect of a larger off-axis fleld which has cosh-like 
dependence in the vertical plane for the planar undulator or in both planes for a helical 
undulator. 

In other words, the most general form of the Wigner distribution obtained from incoherent 
addition of radiation from individual electrons is of the form 

Wir, e) = Nej Wo(r, 0; Ve)P(Ve) rf^Ve, (98) 

Here is the total number of electrons inside a bunch described by the probability density 
function P{Ve) such as J P(Ve)rf^Ve = 1 with Vg representing 4 phase-space variables of 
electrons: two transverse components of position (xe,?/e) and angle (Oex^Oey). If, on the 
other hand, the different electron trajectories simply lead to an offset in position and angle 
of the WDF of a single electron 

Wo{t, 9; re, d^) = Wo{t - Fe, 9 - 9e), (99) 

the integral of Eq. [98] is then replaced with a convolution integral 

Wir, Q) = N, (Iwoir- r,, 9 - ee)P(re, 9,) Sv, SQ,. (100) 



The effect of energy spread in electron beam can be quite signiflcant, and most generally 
it is accounted by extending the integration variable Vg to also include the energy. E.g. the 
effect of a small energy spread 5^ = A7e/7e <C 1, where 7e is the normalized electron energy 
(later denoted as simply 7), leads to 

W{v, Q) = N,jj W^{v - re, 9 - 9e; 5,)P{v^, 9e, 5,) dh^ d'e^ dS,. (101) 

Note that for a small energy change ^ i/Nu with a large number of undulator periods 
Nu ^ the effect on the radiation pattern at a given frequency ujq is identical to that of 
the on-energy particle AS^ = while tuning the radiation frequency off the resonance by 

Auj/uJq — 25e. 



The evaluation of Eq. 101 can be quite involved in terms of computational resources 
required even if being straightforward in all other respects. However, if the electron distri- 
bution P{Ve) is separable, i.e. P{Ve) = Px{xe^Oex)Py{yeiSey)P-f{^e)^ then the 2D projection 
of the WDF, Bx{x^ 9x) and By{y^ 6y) can be easily computed. 
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It is instructive to consider the requirements for when Eq. [99] is apphcable in case of a 
planar undulator (horizontaUy deflecting). As mentioned previously, two eflFects can change 
the shape of the WDF depending on the (small) electron trajectory offsets in position and 
angle in vertical plane. One is the cosh-dependence of the vertical fleld, whereas the other 
is the natural undulator focusing. 

The equation of motion for the average vertical position yav when B^^ = for the undu- 



lator, Eqs. 97, can be written as [33j 

'''^^'^^ = -klVav, (102) 



where the vertical focusing strength is given by k^y = Boye/y^jnieC^ or in terms of the period 
of oscillations due to focusing = 27r//c^^ = V^^Xu/ Ky^ where Ky is the undulator K- 
value and 7 = E/rrieC^ is the normalized energy of the electron. Typically, L^^^ ^ L^^ i.e. the 
slow oscillation phase increment due to the focusing is <C 2ti in undulators. Nevertheless, 
in order to be able to treat vertically oflFset trajectories as simple copies of each other, we 
require that the slow sine-like oscillations due to focusing produce a change in the electron 
trajectory's deviation over the length of the undulator that is much smaller than the natural 



cone of the radiation, yJXjL^ [6j. Integrating Eq. 102 for a typical vertical size Uy^ the 
angle change of the electron trajectory is of the order ayk^^Lu^ which leads to the following 
requirement 

Similarly, the vertical dependence of the magnetic fleld in the undulator By oc cosh(/c^iy), 
leads to the vertical trajectories with an oflFset to effectively sample a larger Ky value. 
Therefore, to enable the simpler treatment, we require that ^KyjKy ^ [kyyY /2 produces 
a change in the undulator wavelength A = A^i/27^(1 + Ky/2) which is much smaller than 
the natural undulator bandwidth A A/ A ^ 1/Nu- This leads to another requirement for the 
electron beam size 



A^ / 1 /2 + X2 

Electrons coming with a vertical angle into a planar undulator generally sample a more com- 
plicated magnetic fleld pattern, such as shown in Fig.[7| Therefore, the following requirement 
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FIG. 7. An example of a vertical trajectory in a planar undulator and the magnetic field as seen 
by the particle. The angular offset is taken to be rather large to illustrate the effect of cosh-like 
dependence of the field on vertical postion. 



can be imposed on the vertical angular size ay 



2 + 



Cfy' <C 



2^K 



3/2 



2 + ^2 



(105) 
(106) 



In summary, if the requirements 103 , 104 , and 106 are satisfied, the simple convolution of 
a single electron radiation pattern with that of the electron bunch phase-space distribution, 
Eq. |100 or Eq. |101 can be used. Otherwise, the more general integral, Eq. |98j needs to be 
evaluated. We note that the potential complications discussed here apply only to the vertical 
plane for a planar undulator with exact translational symmetry of fields in the x-direction. 



C. Revisiting emittance definition 

Rms emittance, Eq.[T4| is widely used as a measure of beam quality in accelerator physics. 
While this definition is attractive due to the fact that it can be applied to a variety of different 
distributions, the connection of the rms emittance to phase space density or brightness 
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available in the beam is generally distribution dependent. Whereas equilibrium processes 
(e.g. radiation damping in storage rings, equilibrium beam in a focusing channel under the 
influence of space charge [34j, etc.) lead to a Gaussian distribution in phase space, beams 
in linear accelerators are rarely in equilibrium. As a result a meaningful characterization of 
the phase space of electron beams or, as we shall see later, the synchrotron radiation, needs 
a more flexible metric than the rms emittance alone. Short of the complete knowledge of the 
actual phase space distribution, a useful way to reduce and represent the information is to 
extend the concept of the rms emittance to the so-called brightness curve or rms emittance 
vs. beam fraction [35j. As we will see, a wide class of practical phase space distributions can 
be effectively characterized by such a curve as the beam fraction is varied from to 100%. 
Three parameters, the usual rms emittance (e = e(100%) with 100% denoting that the 
entire beam is included in the emittance calculation), core emittance, 6c, and core fraction 
fc can convey the information not only about the second moments of the beam distribution, 
but also the peak brightness and what fraction of the beam effectively contributes to this 
brightness. The situation is somewhat analogous to how the peak height and the full- width 
at half maximum complement the rms width information for arbitrary (unimodal and flnite 
integrable in the second moment sense) pulses. 

Below is one prescription for obtaining emittance vs. fraction curve. Here we only consider 
the case of a two-dimensional phase space, x = {x^py where x is the transverse coordinate 
and p can represent (normalized) transverse momentum or angle. The phase space distri- 
bution function P{x^p) is assumed to be normalized, JJ P{x^p) dx dp = 1. One can apply 
the following procedure: 



1) For an ellipse of a flxed area vra, choose Twiss parameters T of the ellipse (c.f. Eq. 14) 
that maximize the beam fraction contained therein: 



/(a) = max / / P{x,p) dxdp with D{a) = {x : xTT"^x < a}. (107) 

iDia) 



2) Obtain the rms emittance e{a) for x G D{a) of Eq. 107 



/ f f P(x v) 

e(^) = V {^^)d {p^)d - {^p)d^ with {u)j^= / / u \ dxdp. (108) 

The parametric curve {f{a)^e{a)) is the emittance vs. fraction curve, e{f). 
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3) Define the core emittance, ec, and the core fraction, /c, according to 

de{f) 



(109) 



df 

fc : e(/c) = ee. (110) 

We have assumed that each individual eUipse D{a) remains centered around the origin as 
does the corresponding centroid of the beam fraction. Generahzation to when this is not the 
case is straightforward by aUowing the chpping elhpse to shift. This procedure for obtaining 
emittance vs. fraction curve is meaningful for distributions which are unimodal (i.e. with a 
single hump) and finite integrable (for second moments). 

It is easy to show that the core emittance is directly related to the peak phase space 
density or brightness Pq = max{P(x,]9)}: 

To see that one simply needs to note that a small area clipping ellipse in the limit a ^ 
cuts out a uniform slice containing the beam fraction naPo and having the rms emittance of 
a/4. It is interesting to note that because of the Property |3| which states that a maximum 
Wigner distribution is /i/2 for any even pure state, and a corresponding 2D equivalent in 
optics of A/2, the minimum core emittance (the diflFraction limit) is therefore 

min(e,) = ^, (112) 

and it can only be larger for a symmetric mode when the coherence /x^ < 1. Thus, the 
core emittance (or peak brightness) is a more general than the rms emittance indicator of 
whether the radiation is coherent. This is because the rms emittance minimum is restricted 
only to a Gaussian coherent mode, whereas the minimum core emittance is realized for any 
symmetric coherent mode. 

Fig. [8] further illustrates the concept of the emittance vs. fraction by showing the curve 
for 3 different distributions: uniform, gaussian, and elliptical. The correlation in x and p is 
removed and the units for x and p are chosen so that the distributions can be written as a 
radial function of r = y^x^~+p^. Furthermore, to facilitate the comparison, each distribution 
is normalized to have e = = ap = 1 in these natural units. As seen, the core emittance 
conveniently captures the fact that the peak brightness of a Gaussian is x 2 larger than that 
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FIG. 8. Radial phase-space distributions (left) and corresponding emittance vs. fraction curves 
(right). All distributions are scaled to have e = 1. Core fraction and emittance for different 
distribution types are shown as well. 



of the uniform distribution of the same rms width, as well as the fact that the core fraction 
in the Gaussian is smaller (0.715 vs. 1 for the uniform). 





e = 0.5 



0.18 



FIG. 9. Radial phase-space distributions (left) and corresponding emittance vs. fraction curves 
(right). Both distributions are scaled to have e = 1. Core fraction and emittance for the two 
distribution are shown as well. 



Another example, Fig. 



9l compares a Gaussian distribution P(r) = 2 and P(r) = 
^^^^-e 2ei + ^^2, with p = 0.5, 61 = |, and 62 = |. The total emittance in this case 

is 6 = pei + {1 — p)e2 = I. Once again, the information about the peak brightness is lost 
with the rms emittance quoted only, but is conveyed conveniently with the three parameters: 
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{e, ec, fc}- A practical measure of beam brightness available can be defined as fc/ ec, a subject 
that we explore further below. 



D. Possible definitions of brightness 

As a phase space quasi-probability, the WDF is the generalized brightness (also known 
as microscopic brightness [35j), i5(r, 9) = W{r^ 9). It is convenient, however, to be able to 
reduce the information to a single parameter, which, for example, can facilitate comparison 
of various partially coherent sources. Here we revisit several of the definitions that can be 
useful for this purpose remembering that no single reduced parameter or a figure of merit 
can suit all the practical purposes. 

1) The following definition, which we denote as classical, can be written (modulo a prefactor 
that generally depends on the actual distribution shape) as 

Bel = (113) 

T is the overall (spectral) fiux. In the definition above we have assumed that the 4D 
emittance can be represented as a product of two 2D emittances. This definition, which 
gives a positive quantity, is easy to compute and can serve as a measure of brightness. One 
drawback is in the use of rms emittance, which, as discussed previously fails to capture 
the peak brightness available in the beam and tends to exaggerate the importance of 
tails when non-Gaussian distributions are encountered. A possible modification to the 
definition of Eq. |113| can be made to write the effective brightness in terms of 



Ba^ait = (114) 

^cx^cy 

where in place of rms emittances as a measure of effective phase space area we use the 
core emittance ecx,y while at the same time reducing the participating fiux by the product 
of the core fractions in each plane fcx,y All the necessary quantities in Eq. 114 can be 
found from the WDF as discussed previously. 

These classical definitions, however, fail to capture the concept of coherence. A mode 
with a large dispersion (emittance) but perfectly coherent is indistinguishable from its 
incoherent analog of the same emittance. 
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2) As we have seen, the WDF contains the information about the density matrix, which, to 
the overaU flux factor leads to the foUowing natural deflnition for brightness (denoted as 
average brightness) 



Bav = ^ J . (115) 

This deflnition is discussed in a classical context in [36], though its genuine justiflcation 
becomes clear from the connection to quantum or wave phenomena. The brightness of 



Eq. |115| is higher for more coherent radiation, even though the dispersion or emittance no 
longer comes into this deflnition. In particular, as pointed out previously, a pure mode, 
no matter how dispersed it gets, would have the same Bav provided the flux remains 
unchanged. 

3) Another deflnition is simply to quote the on-axis peak brightness 



Bo = W{T = 0,e = 0). (116) 

An obvious drawback of this deflnition is that the WDF is not guaranteed to be posi- 
tive. However, as previously discussed, the on-axis WDF is always positive for symmetric 
(even) modes, which are of most practical interest for synchrotron radiation. Addition- 
ally, the peak brightness due to the boundness property (Property |3]) can serve as a 
measure of coherence because for any pure and symmetric (even) mode Bq is guaranteed 
to be related to the total (coherent) spectral flux according to 

Bo,pure,ev ~ \^ ) ^ ' (H'^) 

As pointed out previously, the core emittance is inversely related to the peak brightness. 

Finally, for the purpose of the numerical examples below, it will be convenient to con- 
sider 2D projections of the WDF which are easy to visualize. The extension of the above 
deflnitions to 2D is straightforward and the equivalent meaning remains intact only when 



the mode is separable in x, ^/-planes. In particular, the Eq. |115 in 2D becomes 



B„„ . um^. („8) 
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and equations Eqs. |116| and |117 



^o. = ^x(x = 0,^^ = 0), (119) 
2 

^Ox,pure,ev ^V*^ ' (120) 



with equivalent expressions for ^-plane. We are going to use Eq. |118 even when the mode 



is not separable as a measure of effective average brightness in one plane. 



E. Numerical examples 

Here we demonstrate numerical examples of using the Wigner distribution function for- 
malism for undulator radiation. We start out with a zero emittance electron case. For 
convenience, we scale the results to 100 mA average current for otherwise perfect (pencil) 
electron beam. Throughout all the examples, the electron energy is set to 5GeV and the 
undulator period A^^ = 2 cm. 



Fig. [To] illustrates the calculated angular flux and the central cone for an undulator with 
A^^^ = 250 periods. To convert the computed from radiation flelds quantities to the standard 
units of photons/s/0.1%BW for spectral flux and corresponding angular (per mrad^ or mrad 
in one plane projection) and areal densities (per mm^ or mm in one plane projection), we 
note that the spectral flux density is related to the computed flelds |E(r)p in frequency 
domain according to 

|E(r)p, (121) 



dAduj/uj e nh 

where / is the average beam current (non-FEL process is assumed), and cq is the vacuum 
permittivity. 



Fig. 103. checks the angular flux density (or the on-axis radiation fleld) against the well- 



known expression for the planar undulator 



d'J^ I 



= -aN^YKiK), (122) 



dQdu) juj e 

with flne-structure constant a, and function Fn{K) = K^n^/(1 + ^y[JJ] for the harmonic 
number n and the undulator K where [JJ] = [J{n-i)/2{0 ~ -^(n+i)/2(C)]^ terms of the Bessel 
functions and ^ = nK'^ /{4: + 2K^). The area dA and solid angle dVt elements are related 
via dA = R^dVt with R being the distance to the source, allowing to cross-check the angular 
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FIG. 10. Comparison of calculated (dots) vs. theoretical values (solid curves) for undulator ra- 
diation: (a) angular flux density, (b) central cone flux, (c) on-axis 4D brightness, (d) on-axis 2D 
brightness. The first three odd harmonics are shown. The undulator is planar with Nu = 250. 
Refer to text for other parameters. 



flux, Eq. |122[ in terms of the computed fields via Eq. |121[ Fig. [I0| 3 compares on-resonance 
spectral flux with the analytical result 



J- II 

-—^ = --TiaNuQn{K), 
d(jj/(jj 2 e 



(123) 



where Qn{K) = (1 + K'^ /2)Fn{K) /n. In what follows, we denote the spectral flux by 
simply J^o implicitly assuming the usual 0.1% bandwidth scaling. To find the total flux, the 
detector in simulations is placed 50 m away from the undulator center and the electric field 
is computed on a 1024 x 1024 3 mm square grid. 

To check that the code correctly computes on-axis (peak) brightness for 4D and 2D WDF 
computed from the fields, we use the Eq. |123| with Eqs. |117| and |120[ which relate the total 
flux to the peak brightness of any symmetric coherent mode according to Bq = (2/A)^J^o 
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and Box = (2/A)J^o- The results of this cross-check are shown in Fig. 



10 



and 



lOi. 



1. Example: helical undulator on resonance 



Wigner-Stokes distribution function (pli/s/0.1%BW/mm/mrad) 
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FIG. 11. Wigner-Stokes density distribution functions computed for a helical undulator. The x-ray 
phase space is back-propagated to the undulator center. Refer to text for details. 



While the planar undulator radiation on- axis is fully horizontally polarized, Fig. 11 shows 



the WDF for a helical undulator at its first harmonic (A^^^ = 250, Kx = Ky = 0.696, huj = 
8keV). The WDF is obtained from the detector plane placed 50 m away from the undulator 
center, and subsequent back-propagation of the radiation phase space back to the center of 
the undulator. As discussed previously, the case of a (nearly) pure circularly polarized wave 
leads to I^Sol = I^Sil with other generalized Stokes parameters being approximately zero as 
seen in Fig. [TTj 

The x-like shape of the x-ray phase space is persistent throughout all the examples. 
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The explanation behind it is simple — undulator, being an extended source, has radiation 
emitted from its beginning and the end, which must advance different distances to reach 
the observer, or when (back) propagated to the undulator center. This results in the x-like 
shape, with the two branches corresponding to the undulator ends. 

In the remainder of this section, we limit our numerical examples to planar undulators 
investigating x-ray phase space for radiation on and off resonance, the segmented undula- 
tor with a quadrupole focusing in between, and a 25-m long undulator including electron 
emittance and energy spread effects. 



2. Example: planar undulator on resonance 




FIG. 12. (a) Wigner 2D distribution for a planar undulator with — 250 on resonance of the 
1st harmonic at 8keV photon energy along with (b) the emittance of light vs. fraction curve. The 
WDF is back-propagated to the undulator center. 



This example illustrates the x-ray phase space for radiation at undulator resonance. 
Fig. 12, along with the emittance vs. fraction curve. It is seen that > 1 or emittance 



is not the minimum possible for the fully coherent mode. On the other hand, the core 
emittance is its possible minimum as discussed previously. Also, note the value of the /3- 
function or Rayleigh range, is somewhat different than 1/^^/2 or I/^/27r values commonly 
quoted in the literature. Additionally, the full beam and its core have different /3-function 
values. Therefore, a proper matching with the electron beam depends on whether one 
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maximizes the peak brightness or minimizes the overaU rms emittance of hght. 



3. Example: segmented undulator with quad focusing 



Next, we consider a segmented undulator with a quadrupole focusing in between the two 



segments. Fig. [13] shows trajectories for different horizontal offsets of electrons going into 
the undulator. 
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FIG. 13. Trajectories inside a segmented undulator with a quadrupole focusing. The undulator has 
two segments each with — 100 periods. The separation between the two is 0.486 m. A single 
horizontally focusing quadrupole of length 0.3 m is located at the center with 3.5 T/m gradient. 



Fig. 14 shows the flux at a detector 30 m away from the center of the undulator and 
the WDF back-propagated to the undulator center. It is seen that in this case the electron 
interferes with itself and the WDF clearly shows features present in a coherent superposition 
of two modes. 
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FIG. 14. Radiation of the 1st harmonic on 8keV energy resonance from the segmented undulator 



with quadrupole focusing of Fig. [T3j Different rows correspond to different trajectory offsets, Xini, 
as shown. The left column is 2D WDF back-propagated to the center of the undulator, which 
corresponds to the right column showing the radiation spectral flux density at the detector 30 m 
away from the undulator center. 



4. Example: radiation off undulator resonance 

Here we consider the radiation off the undulator resonance. This is not only of interest 
for practical cases of detuning or selecting photon energy in a monochromator but also when 
considering off-energy electrons (electron beams with energy spread). This is because for 
undulators with a large number of periods, the effect of tuning off resonance is identical to 
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keeping the radiation frequency Uq the same but changing the electron energy according to 
Auj/ujq — 2A7/7. 




N^(Aco/cOq) 



FIG. 15. Scanning around the 1st harmonic resonance of 8keV for A^-^^ = 250 period planar 
undulator: (a) angular flux on axis, (b) integrated (central cone) flux. Red circles denote analytical 
values. 



Fig. 15 shows the effect of the radiation frequency scanning around the resonance of the 
1st harmonic on the angular (spectral) flux density on axis and integrated flux. Twice the 
spectral flux is available for radiation Auj/ujq ^ 1/Nu below the resonance. 



Fig. [16] demonstrates the light phase space for 3 different values of the radiation frequency 
detuning along with the radiation pattern (50m from the undulator center). 



Fig. 17 shows the omittance and /3-function of light for scanning the radiation frequency 
around the resonance. As shown previously, the core omittance is A/Stt in all cases, whereas 
the rms omittance is minimal (though with > 1) around the resonance. 



Finally, Fig. 18 shows the effective 2D average brightness Bavx = JJ W^{x^9x) dx dOr^ 



and fi^g^ = ^ If ^xi^^ ^x) dxdOx where the normalized WDF Wx{x, Ox) = Bx{x, 9x)/J~' with 
= JJ Bx{x^9x) dxdOx- The deviation of /x^^ from 1 is due to the fact that the radiation 
mode is not separable, even though the radiation is fully transversely coherent in this case 
and therefore the full 4D = 1. The peak 2D brightness, which is not shown, simply 



follows the trend of Fig. 15 since it is related to the flux according to Bqx = T{2/\). 
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FIG. 16. The WDF (top row) corresponding to Fig. 15 for different detuning off the resonance. The 
WDF is back-propagated to the undulator center. Sohd hnes show position and angular intensity 
projections. The radiation pattern used to calculate WDF (bottom row) is obtained at 50 m from 
the undulator center. 



5. Example: including emittance and energy spread of electrons 

Here we provide an example of including emittance and energy spread to the calculated 
WDF. For simplicity, we continue to limit ourselves to 2D projection of the Wigner distribu- 
tion function and treat electron phase space probability distribution function as separable 



P(re,0e,^^e) = Pxi^e^ ^ex) PyiVe^ ^ey) Ps{^e) - Fig. 19 shows the horizontal phase space at 
5GeV obtained from the simulations of the photoinjector for 77 pC per bunch and 1.3 GHz 
repetition rate (average current of 100 mA), including the effects of the merger and the linear 



accelerator [37j. See Fig. 19, The energy spread of the electron beam is as^ = 2 x 10~^. To 
illustrate its effect, we consider a 25-m long undulator with A^^^ = 1250 periods. Table [TT] 
summarizes the parameters used in this example. As seen, the radiation is computed slightly 
below the resonance where the flux is roughly doubled. 

To provide more optimal matching for the core of the beam, f3x is chosen to be /3a; = 4m 
close to the Rayleigh range of the core of the radiation from a pencil (zero emittance) beam, 
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FIG. 17. Light 100% and core emittances (left) and beta function (right) at the undulator center 



corresponding to Fig. 15 for different detuning off the resonance. 




Fig. 20=1. Fig. 20 3 shows the effect of the energy spread for otherwise ideal (zero emittance) 
beam. Some degradation of the Bavx can be seen. 



Fig. [21] shows the effect of the beam emittance on the radiation phase space. It can be 
seen that the negative parts of the WDF are no longer present after the convolution, and 
the average 2D brightness is degraded by a factor of about 2 and 3 for 25 and 100 mA cases 
respectively compared to zero emittance. Other relevant parameters of the radiation are 
shown in Fig. I2l} 
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FIG. 19. Electron phase space at the center of 
currents: 25mA (left) and 100mA (right). 




m long undulator corresponding to two different 



TABLE II. Parameters used in computing the radiation phase space. 

Number of periods, Nu = 1250 

Undulator period, A^^ = 2 cm 

Harmonic number, n = 1 

Resonant photon energy, hjj — ^ keV 

Detuning radiation frequency, Acj = —{).7^uj/Nu 

Beam energy, E — ^ GeV 

Electron energy spread, a^^ = 2 x 10~^ 

Electron emittance, ex — 11,29 pm 

Average current, / = 25, 100 mA 

/3-function, /3a; = 4m 



V. CONCLUSIONS 



The Wigner distribution function approach to describe partially coherent radiation in 
phase space has been presented. Despite the general nature of the approach, the true power 
of the method to describe modern and future x-ray synchrotron sources is to employ 5 or 6D 
phase space (times 4 for arbitrary polarized light) complementing the 4D transverse phase 
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FIG. 20. X-ray phase space corresponding to zero emittance and zero energy spread (left) and 
non-zero energy spread as^ = 2 x 10~^ (right). The beam current is 100mA. 




FIG. 21. X-ray phase space including the effects of emittance and energy spread for the two 
different currents. 



space with frequency and time where the timing structure is important. Though straightfor- 
ward, such a description is rather challenging from the point of computational requirements, 
even though a sampled approach similar to particle tracking in accelerator physics can be 
employed to represent the radiation in the entire 6D phase space (the microscopic brightness 
is allowed to take on negative values). When the x-ray optics beamline consists of drifts and 
perfect lenses without clipping apertues, this description is complete and allows to fully ac- 
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count for the light properties foUowing geometric optics transformation rules. Introduction 
of apertures in the beam, however, requires the convolution of the transmissive mask's WDF 
with that of the beam. In this cases, it might be more efficient to consider decomposition of 
the partially coherent light into orthogonal mutually incoherent modes and to include the 
diffraction effects on each mode separately. 

Nevertheless, the Wigner distribution function is demonstrated in this paper to be a 
rigorous and insightful way to describe the coherence and other properties of the synchrotron 
radiation. Its use will grow in importance as synchrotron x-ray sources with higher coherence 
become more prevalent. 
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